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ABSTRACT 


A waterflood of the Second Wall Creek Sand of Naval Petroleum 
Reserve No. 3 is being considered in order to increase the ultimate oil 
production. A streamline model of the Second Wall Creek Sand was 
developed here and used to predict the response of the reservoir toa 
waterflood. 

Results of the model study indicate that the oil production of the 
Second Wall Creek Sand can be increased by as much as 15, 000, 000 
barrels over a twelve year period by waterflooding. This would bring 
the total production available from the Second Wall Creek Sand to around 


32,000, 000 barrels. 
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CHAPTER I 


INTRODUCTION 


Naval Petroleum Reserve No. 3 is located in Natrona County, 
about thirty-five miles north of Casper, Wyoming. The major geo- 
logical structure on the reserve is an anticline called ''Teapot Dome". 
The Second Wall Creek Sand is the most important of the several oil 
bearing formations on the anticline. 

The Second Wall Creek Sand is a solution gas drive reservoir. 
There appears to be very little natural water drive. Early estimates, 
based on recovering one-seventh of the original oil in place, predicted 
the recoverable reserves from this formation to be around 17,000, 000 
barrels(1). 

Waterfloods of the same formation in neighboring fields have 
proven to be successful in increasing production. Therefore, in order 
to enhance the recovery, the development and production of the Second 
Wall Creek Sand as a waterflood is being considered. 

The objective of this work is to apply streamline simulator 
techniques to the Second Wall Creek Sand and predictthe results ofa 


waterflood under certain operating conditions. 





CHAPTER II 


THE STREAMLINE MODEL 


The theory and development of the streamline model have been 
rigorously described by LeBlanc and Caudle(2), Rust(3), Wessels (4) 
and Kazmann(5). A review of the assumptions and the mathematics 


Psepresentea here, 


oe Basic Assumptions 
The basic assumptions used in formulating the streamline model 
are as follows: 
1, Negligible gravitational effects within the reservoir. There 
is no gravity segregation of fluids. 
2. Steady-state flow within the reservoir. Incompressible 
reservoir fluids. Operations are at nearly constant pressure. 
3. Homogeneous, isotropic reservoir. 
4. Uniform dispersion of free gas. 


5. Onilbank formation. 


The simulator used is designed specifically for 011 bank build-up 
and will not accurately predict the results of other operations such as 


pressure maintenance. 





ae Fundamental Equations 


The starting point for the development of the streamline model is 
the continuity equation. This partial differential equation is derived from 
a three dimensional material balance and is written as: 


af 3 3 


The volumetric flux may be expressed by Darcy's Law in the 
Sitterential form as: 


k d® 
= - = (2.2.2) 


eee o 
A A 
When the volumetric flux components in the x, y and z directions, as 
described by Darcy's Law, Equation 2.2.2, are substituted into 


Equation 2.2.1, the following equation results: 


Oe ae 3 a 


3 ( . EEN 30 
Ox Ox oy 


0M; Oz WU 382 ot (2-255) 
By using the assumptions listed previously, Equation 2.2.3 
may be simplified. The resulting equation may be written in the form 
of a La Place equation, 
Z 
ap oP 


oP a 4 = @) (2.2.4) 





The solution of Equation 2.2.4 yields the pressure at any point 
in a porous medium for which the assumptions made are valid. Knowing 
the pressure at all points in a porous medium, the average velocity ofa 
fluid particle, and thus its movement, may be computed. 

By transforming to cylindrical coordinates and allowing radial 
flow about a point, Equation 2.2.4 may be solved for the flow potential 
in a horizontal, isotropic, homogeneous porous medium. Restoring 
the resulting equation to rectangular coordinates and using the rule of 
superposition for a multi-well system, the solution of Equation 2.2.4 


becomes 


poor) = my - ge Dam [Ox ey) | e290 





where p_ is the mean reservoir pressure. 
m 


In a porous medium, the relationship of flux to velocity is 
dp 
hell al ITS: 


By differentiating Equation 2.2.5 with respect to x and substituting the 
resulting derivative into Equation 2.2.6 the velocity in the x direction 
in a rectangular coordinate system becomes: 


> ale 


i n 
v. &, y) = - inh = qd. Gay ae (cuc. t} 





In a similar manner, the velocity in the y direction becomes: 


l = YOY; : 
a > | 2 | ier oe eS 

eu a ae 
y th i an 5 (y-y,) (2.2.8) 


Using vector addition, the velocity of a particle of fluid is: 


oe ae (2.2.9) 


Zn. Finite Difference Equations for Movement 

In order to apply Equations 2.2.7 and 2.2.8 to the streamline 
model, it is assumed that the velocity remains constant during any one 
small distance increment, 4s. Then the time increment, At, is 
calculated by: 


A 
ae ae (cca) 
V 


From the particle's starting point (x , y ), it will move to a new point 
s Ss 


calculated as follows: 


= + F A 
eee #5 ee - y,) : 


(Ze 3 eZ) 


Ys+1 ime ‘y Se Y slat 

By starting on the wellbore of an injection well and successively 
applying Equations 2.2.7 through 2.3.2, a trace of any particle's move- 
ment may be made. This trace, called a streamline, is followed until it 


reaches a production well whereupon it is terminated. The travel time 


fora single streamline is the summation of the time increments 





corresponding to the number of steps taken along the streamline. 


2.4 Streamtube Concept 

In actual computation, a representative number of streamlines are 
chosen to represent all fluids emanating from the source. An example of 
the streamlines for a bounded pattern are shown in Figure l. A few 
typical streamlines are shown in Figure 2. By dividing the distance 
between two adjacent streamlines into two equal parts, the streamline 
may be said to represent the flow of all the particles within the dividing 
lines. A streamtube has thus been described and is shown in Figure 3(3). 
As a fluid particle is traced along its streamline, the volume of fluid 
produced at the streamline's production well can be calculated and 


accumulated. 


25.5 Multiple Fluid Flow 

The streamline model described so far deals with a single fluid. 
In a waterflood operation, water is injected into the reservoir to displace 
the oil and gas. Therefore as many as three fluid regions may be present 
in the reservoir at any time and may each have a different mobility. 

The first assumption for the multiple fluid system is that the oil 
is miscibly displaced by the water in a piston-like manner. The oil is 
banked up ahead of the water such that all oil displacement takes place at 
the leading edge of the water front. A diagram of the concentrations is 


shown On page 9. 








@® Injection Well 
O Production Well 


Figure 1 Streamlines for a Bounded Pattern. 
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Figure 2 Typical Streamlines. 
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Figure 3 Typical Streamtube. 
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Water ed-out Pre-flood 
Region Conditions 


With the assumption of piston-like displacement the displacing fluid inter- 
face can be easily followed along the center pathline of a streamtube. 

The most significant approximation used in extending the stream- 
line model for multiple fluids is that the pathlines will be the same as the 
streamlines generated for a single fluid. This approximation, while not 
correct, has been found to be reasonable except at extreme values of 
mobility ratio(3). 

As the fluids proceed along the streamtube the flow rate of the 
fluid stream will change because of the different fluid region mobilities. 
To correct for the changing fluid flow rate the concept of conductivity 
ratio is applied. Conductivity and conductivity ratio are developed and 
explained by Caudle(6). Briefly the conductivity ratio is a comparison 
of the conductivity of the medium at some time when more than one fluid 
is flowing to the conductivity of the medium initially when only one fluid 


was present. The general equation for the conductivity ratio is 





10 


y =— (2.5.1) 


k=1 : k 
The velocity Equations 2.2.7 and 2.2.8 are multiplied by the conduc- 


tivity ratio to correct for multiphase flow. 


Zo Modification for Thickness Variation 

An empirical modification to allow for varying thickness has been 
developed by Wessels(4). This modification is useful if the thickness at 
each point in the reservoir can be defined. A system of grid squares 
with assigned thicknesses is used to describe the reservoir thickness 
variation. 


The modifications for thickness variation are as follows: 


1 > q. x -X, on 
ae 6 h. 2 Pale he ee 
i=l 1 (sx -x,) + (y-y, ) j 
1 i q. y-Y; Ds amie: 
as —-_ 7 u-_ ° r 
Re fs |e GaP 


The factor tne is an approximation factor which reflects an increasing 
J 

velocity as the thickness decreases and a decreasing velocity as the 

thickness increases. The diagram below shows the relationship between 


reservoir thickness and velocity for a single streamline. 





il 


Production Well 


Injection Well 





Front Location (x, y) 
Thickness nh; 


Paes Mathematical Bounding of the Streamline Model 

An imaging technique developed by Lin (7) which can determine 
the necessary rates for a preselected set of image wells that will 
confine a system of streamlines within a given area was used for this 
study. The technique uses a least squares method to determine the 
image well flow rates so that a no-flow boundary is established between 


each of the reservoir bounding points. 





CHAPTER III 
APPLICATION OF THE STREAMLINE MODEL TO THE 
SECOND WALL CREEK SAND 
aa) Reservoir History 

President Wilson's Executive Order of April 30, 1915 designated 
the Teapot Dome area in Wyoming as Naval Petroleum Reserve No. 3 
(NPR 3). On April 7, 1922, the reserve was leased to the Mammoth Oil 
Company for the purpose of exploitation. Following a decision of the U.S. 
Supreme Court all of the producing wells on the reserve were shut-in on 
December 31, 1927. The field remained shut-in until 1951-1953 when an 
exploratory program was initiated; thereafter it was shut-in until 1958 
when an offset drilling program was instituted to protect against drainage 
by adjacent operators. 

The wells are distributed unevenly in the Second Wall Creek Sand 
because development was started at the north end of the reserve and was 
stopped prior to completion. Based on approximate surface area the 
northern third of the reserve has 69 wells, the middle third has 23 wells 
and the southern third has 13 wells. 

Approximately 3, 600,000 barrels of oil were produced prior to 
December 31, 1927. 1,200,000 barrels of 011 have been produced since 


1951 with most of this production occuring after 1958. 
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Be: 


Bec Location of Oil and Gas 

Figure 4 shows the approximate locations of the gas cap, oil zone 
and water in the Second Wall Creek Sand (8). 

The Second Wall Creek Sand just beyond the northwestern boundary 
of NPR 3 has been under waterflood for several years. The Navy has been 
producing wells along this boundary to protect against drainage. Conse- 
quently a portion of the Naval Reserve has been swept by the waterflood. 
To allow for this production, 170 acres along the northwestern boundary 
were omitted from the area which was modeled. The total area modeled 


was 3000 acres. 


oie, Reservoir Properties 
Reservoir data and operating conditions were obtained from 
Tesoro Petroleum Corporation and U. S. Navy drawings of NPR 3 
(8, 10, 11), Geological Survey and U. S. Navy publications (1, 9), 
U. S. Navy correspondence (12,13) and various core data and logs pro- 
vided by the Officer in Charge of Naval Petroleum Reserve No. 3. 
Table I lists the reservoir and fluid properties used for the 
modeling study. A discussion of significant individual items follows. 
Porosity and Saturations The average porosity and connate water 
saturation used are the same as those used in previous estimates of 
waterflooding at the north end of the reservoir (12, 13). No significant 


difference was noted between these figures and the core data examined. 
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Figure 4 Location of Oil and Gas in the Second Wall Creek Sand. 
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RESERVOIR AND FLUID PROPERTIES 
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For this study it was assumed that the waterflood would start at 
the very beginning of the production. Whe model was run with an initial 
gas saturation of 5% to allow oil bank build-up to occur and to show the 
effect of fill-up time. The initial oil saturation used was 60%. 

The residual gas saturation was set at zero because it was 
assumed that all of the free gas would either be displaced or put back 
into solution as the oil bank passed. 

The residual oil saturation after the waterflood passed was set 
at 20%. This is almost twice the 12% average residual oil saturation 
found in the cores available. The residual oil saturation in the cores 
was considered an absolute minimum which would be difficult to achieve 
Putne reservoir. 

Fault Blocks The Second Wall Creek Sand is heavily faulted and 
was, for the presentation of well productivity data, divided into about 
twenty possible fault blocks(12). In order to model the field, taking into 
account the faulting, the field was divided into ten sections along the . 
largest faults. It was assumed that the faults were sealed and the 
sections were then modeled independently. This had the added advantage 
of making the modeled reservoirs smaller and easier to manage. Also, 
properties which vary within the reservoir could be easily changed from 
section to section. Figure 5 shows the configuration of the ten fault 


blocks chosen. 
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Figure 5 Fault Blocks Used in Modeling the Second Wall Creek Sand. 
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Net Sand Thickness Thin beds of shale are common throughout 
the Second Wall Creek Sand, but they thicken toward the south and east 
and split the sand into several distinct layers. Some wells have three or 
more layers of sand separated by shale beds from two to twenty-five 
feet thick(12). 

Combining the information available from the cores with the 
general trends mentioned above, a percentage discount factor was used 
for each fault block to reduce the gross sand thickness toa net sand 
thickness. Table II shows the discount factors that were applied. 
Figure 6 shows the gross thickness of the Second Wall Creek Sand(10). 

Permeability The temperature, pressure and permeability 
were used to determine the fluid and fluid flow properties. The perme- 
ability itself is not used in the streamline model calculations because the 
calculations are made at a constant flow rate. This means that the 
injection pressure will have to be great enough to make the injectors 
flow at the required rate. Rate selection will be discussed in the 
next Chapter, 

Gas-Oil Ratios Gas-oil ratios vary greatly depending on where 
in the field the well is located. The gas-oil ratio and solution gas -oil 
ratio are used by the computer program to indicate the arrival of the 
oil bank. Representative values of the two ratios were used but they do 


not represent a rigid material balance. 
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TABLE II 


DISCOUNT FACTORS APPLIED TO 
GROSS SAND THICKNESS 
TO OBTAIN NET SAND 
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Figure 6 Isopachous Map of the Second Wall Creek Sand- 
Gross Sand Thickness. 
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3.4 Operating Conditions 
Table III is a listing of the operating conditions used for the 


simulation. 


TABLE Itt 


OPERATING CONDITIONS FOR THE WATERFLOODING 
OF THE SECOND WALL CREEK SAND 


No. of Production Wells 35 
No. of Injection Wells 132 
Well Pattern 5-Spot/Line Drive 
Well Spacing, Acres 10 


Maximum Oil Production Rate, 
STB/Well/Day 46-300 


Oil Production Rates During Fill Up, 
STB/Well/Day 


Total Water Injection Rate, STB/Day 





Field Development Based on a proposed Navy Development 
Plan(1l1l) the field was assumed to be drilled on ten acre spacing. 

Well water injectivity was assumed to be in the same range as 
the well's initial oil productivity. Therefore, approximately a one-to- 
one ratio of injectors to producers was maintained. 

The injection-production well pattern used was a combination 


five-spot and line drive. The injection wells were generally placed 
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along the outside, nearest the oil/water contact, in a line drive pattern to 
force the oil toward the center. Then, where area allowed, injection wells 
were placed within the producers in a five-spot pattern. 

Production Well Rates In the streamline model the oil production 
rate of a well during the fill-up time is set equal to the present rate to 
which the well's oil production has declined. After oil bank breakthrough 
the oil production rate of a well is set equal to the initial or maximum 
oil rate of that well. 

The initial rates for each production well were determined from 
a composite drawing of the Second Wall Creek Sand production(8). Since 
the reservoir will be waterflooded from the start, the oil rates during 
fill-up were set at 10% of each well's initial rate. At reservoir conditions, 
the total water injection rate was set equal to the total oil production rate 
after fill-up to satisfy the steady state assumption. 

The individual injection rates were assigned to the injection wells 
in such a manner to contain the streamlines as much as possible to the 
oil bearing area. Some flow into the gas cap and the watertable cannot 
be stopped but proper balancing minimized the losses. 

Figures 7 and 8 show the well pattern used and the streamlines 
generated by the model reservoir. See Appendix A for a listing of the 


individual well data used. 
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Figure 7 Model Generated Streamlines-Northern Half of 
Second Wall Creek Sand. 
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Section 34 
Fault 


X Injection Wells 
O Production Wells 
Scale :1''= 2300' 


Figure 8 Model Generated Streamlines-Southern Half of 
Second Wall Creek Sand. 





CHAPTER IV 


RESULTS OF THE STREAMLINE MODEL SIMULATION 


4.1 Predicted Production History 

Figures 9, 10 and 11 show the predicted response of the Second 
Wall Creek Sand to a waterflood operated as previously described. It is 
predicted that 18,000,000 barrels of oil will be recovered in Fe years 
and 30,000,000 barrels of oil after ten years. Water breakthrough occurs 
around 1,000 days and the water-oil ratio climbs to over five after ten 
years. Ultimate recovery at a water-oil ratio of twenty is predicted 


to be around 32,000,000 barrels and takes twelve years. 


4.2 Effects of Major Assumptions on Recovery Predictions 
Although the oil recovery predicted agrees with the present 
reserve estimate of 28, 800, 000 barrels(9), the predictions of both 
quantity and time are only as good as the assumptions made. 
In order to provide a proper perspective, the variables and 


assumptions which could most effect the results will be discussed. 


Heterogeneity The model was run assuming that the reservoir 
sand in each fault block was homogeneous. Frequently, however, hori- 
zontally bedded sandstones exhibit different permeabilities among the 
individual stratum. The presence of high permeability strata will 


reduce the efficiency of the waterflood. The high permeability zones 
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water-out faster than the rest of the reservoir and then increase the 
water oil ratio much sooner than expected. If high permeability strata. 
are shown to be a significant factor the predictions presented here will 


be optimistic. 


Initial Gas and Oil Saturation The presence of free gas in the 
oil zone of the reservoir means that there will be a certain time lag, 
fill-up time, before the oil bank will reach the production wells. This 
time lag is approximately equal to the time required to inject enough 
water to fill the gas space and the space of the oil that is produced. The 
initial gas saturation chosen to represent the reservoir was 5%. If the 
quantity of gas present in the reservoir is actually less than 5% the oil 
production rate will increase earlier and if the space is assumed to 
be filled with oil, more oil will be recovered. If the gas saturation is 


greater than 5% the opposite will be true. 


Production and Injection Rates The effect of increasing or 
decreasing the injection and production rates would be to condense or 
expand the time scale. If for example it was found that the reservoir 
could operate at injection and production rates 50% greater than those 
assumed here, the recovery that would have taken l, 500 days could be 
achieved in 1,000 days. See Appendix B for an example of rate adjust- 


ment to obtain uniform water breakthrough times. 





CHAPTER V 


CONCLUSIONS AND RECOMMENDATIONS 


See Conclusions 

The waterflood of the Second Wall Creek Sand will increase the 
total oil production to around 32,000,000 barrels. This is 15, 000, 000 
barrels over that which could be expected from solution gas drive primary 
production. At the operating rates used the waterflood will require around 


twelve years. 


Sec Recommendations 

1, It is recommended that the waterflooding of the Second Wall 
Creek Sand be planned and cost evaluations made. The streamline model 
developed here is well suited for studying development and operating 
considerations so that project economics can be evaluated. 

2. Since the effects of high permeability strata on secondary 
recovery operations are always adverse and sometimes severe, it is 
recommended that the possible existence of such strata be thoroughly 


investigated during waterflood planning and development. 
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Symbol 


X,Y,4% 


NOMENCLATURE 


Definition 
Cross Sectional Area 
Reservoir thickness 
Permeability 
Length; number of fluids 
Number of wells 
Pressure 
Flow rate 
Distance increment along path line 
Saturation 
Time step along path line 
Time 
Flux 
Average fluid particle velocity 


Rectangular coordinates 


Conductivity ratio 
Fluid mobility 
Viscosity 

Density 

Porosity 


Potential 
31 





Subscripts 
cw 
gl 


gr 


Ol 


or 


Sif 


Connate water 

Initial gas 

Residual gas 

Index of wells; in 

Index of particle location 
Index of fluids 

Mean 

Sil] out 

Initial oil 

Residual oil 

Starting point of fluid particle 
Production 

Water 

In the x-direction 

In the y-direction 


In the z-direction 
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APPENDIX A 


INDIVIDUAL WELL DATA 
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TABLE IV 
WELL DATA—FAULT BLOCK 1 
iNeter LOumneure TZ 


Well Rate,STB/Day Thickness Well Rate,STB/Day Thickness 





No. (+Inj, -Prod) Fe: No. (+Inj, -Prod) Ft. 
] oie 29 Zu 96 32 
az 54. Ce a2 -~96 a2 
3 -90. 30 23 48 24 
4 SVE 26 24 -72 24 
5 -87. 29 co -72 24 
6 218) Sy 26 -78 26 
ct 56 28 at 84 28 
8 -90 30 28 -90 30 
2 -~96 32 29 96 32 
10 78 26 30 -96 32 
1] -81 ag) Sit -96 a2 
Vz 90 30 a2 48 24 
i -96 32 i) -72 24 
14 96 32 34 72 24 
iS 48 24 35 -72 24 
16 ez 24 36 78 26 
yg -78 26 Shy -84 28 
18 -81 27 38 90 30 
EG 87 a By) -96 a2 


20 -93 al 40 96 pe 
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TABLE IV (Continued) 


Well Rate,STB/Day Thickness Well Rate,STB/Day Thickness 





No. (+Inj, -Prod) Et 4 No. (+Inj, -Prod) rte 
4] 48 24 61 -90 20 
42 -72 24 62 81 cag 
43 ia 24 63 48 24 
44 -72 24 64 36 24 
45 ( 24 65 72 24 
46 -78 Zo 66 -78 26 
47 84 ZS 67 93 31 
48 -90 30 68 -90 30 
49 96 az 69 84. 28 
50 -96 oe 70 -84 28 
Sy 96 32 “1: 48. 24 
Sy4 48 24 2 48. 24 
53 -72 24 chs, -78 26 
54 -72 24 74 -78 Ze 
BS -72 24 75 36 24 
56 90 30 76 40 24 
Si -90 30 1a 48 24 
58 90 30 

Sy) -90 30 

60 90 oi0, 


Production well oil rates during fill-up equal 10 STB/Day 
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FAULT BLOCK 1 


Section 20 
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x Injection Well 


O Production Well 


FAULT BLOCK 2 





Section 29 


Figure 12 Well Locations for Fault Blocks 1 and 2. 
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eile oy 
WELL DATA—FAULT BLOCK 2 
Rerer te higure. [2 


Well Rate,STB/Day Thickness Well Rate,STB/Day Thickness 


Now (+Inj, -Prod) JDie, No. aia, -P rod) Ft. 
1 a 28 20 a 29.4 
Vs “8 3058 21 2 3055 
3 : EM ies, Va fe -175 29.4 
- * Se Za -180 3051 
5 -175 Loe 24 =: ZO 
6 -180 307 25 a Ze 
qd -~195 Blac 26 -155 (da 8, 
8 =2.0 5 34.3 vag -175 (ae oat | 
7) -175 25.6 28 “6 24.5 
1G -185 30.8 ZY -135 22.4 
11 -~200 S220 30 -147 24.5 
cz =2 10 a5 ul 7 Zone 
13 * 29.4 a2 -160 26.0 
14 * 30.8 bie) zs Bose 
15 7 5386 34 -~175 29.4 
16 “ 55% 55 *e Zi 
17 -~175 Zo. 36 ** Z3et 
18 -185 30.8 of, a 2 Sond 
149) -200 Sie 38 * Zo 


Production well oil rates during fill-up equal 20 STB/Day. 


*Program calculates injection well rates based on reservoir thickness 
and steady state material balance. 





fieUs Ira 
WELL DATA—FAULT BLOCK 3 
Metemtorieure 12 


Well Rate,STB/Day Thickness 


_No. (+Inj, -Prod) Ft. 
] 140 ler 
Z 140 35 
3 140 55 
i ~178.5 Sia 
5 -182 36.4 
6 ~178.5 Boe 
7 218.4 36.4 
8 ZO. 2 iit 
2 Z26. 2 Saat 

10 -182 36.4 
1] -192.5 38.25 
V2 -203. 40.6 
13 -188.5 oe at 
14 NT es B5% 

15 182 36.4 
16 203 40.6 
EY -210 42. 

18 -182 36.4 


Production well oil rates during fill-up equal 20 STB/Day. 
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FAULT BLOCK 3 


/ s Ou Fe | : 
x Oo o;1x oOo x 

15 73 il 8S 5 2? 

xX Oy ‘lo x @ xX \ 

16°08 6 
\ x Oy \ Ss Z roy ., 
: \ : . 
Section 33 34 


FAULT BLOCK 4 





Section 33 


FAULT BLOCK 5 


Section 33 





‘ x Injection Well 
O Production Well 


Figure 13 Well Locations for Fault Blocks 3-0. 






TABLE VII 
WELL DATA—FAULT BLOCK 4 


Reter to Pisuneel 3 


Well Rate, STB/Day Thickness 


_No. (+Inj, -Prod) 19 e- 
l 88 44 
Z -100 50 
3 90 45 
+ -100 50 
5 86 43 
6 -94 47 
7 48 48 
8 -~48 ao 
) -48 48 

10 48 48 


Production well oil rates during fill-up equal 10% of initial rate. 
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TABLE VIII 
WELL DATA~—FAULT BLOCK 5 
Refer to ero e V5 


Well Rate ,ote, may Thickness 


_No, (finj, xed) ace 
il -50 50 
Z 50 
3 -50 50 
ot 48 
5 -50 50 
6 * 50 
7 -46 46 
8 46 


Production well oil rates during fill-up equal 3 STB/Day. 


* Program calculates injection well rates based on reservoir thickness 
and steady state material balance. 





Ae ies 
WELL DATA—FAULT BLOCK 6 
Refer to Figure 14 


Well Rate,STB/Day Thickness 


_No. (+inj, -Pred) 1D 
l 260 56 
2 240 Sy e 
3 245 ao 
4 245 53 
5 225 49 
6 -340 5S 
i -300 60 
8 -300 60 
9 -300 60 

10 -285 ii 
ih -2155 Bel 
r2 ou 55 
13 300 60 
14 -300 60 
15 -270 54 
16 ZS 48 
Le 305 5/2 
18 -275 55 


Production well oil rates during fill-up equal 25 STB/Day. 
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FAULT BLOCK 6 


Section 3 : IG 


FAULT BLOCK 7 


he BU Ne ae US 
\ oi 
\ ig 13 \ uf a, \ 
Section 4 vig ig’ \3 x 
8 





FAULT BLOCK 8 


V6 gl [Sale i 
/ Zt It 3 | 
a" o x 
ie 9 / ™~ | | 
a o | 
Section 10 11 


x Injection Well 
O Production Well 


Figure 14 Well Locations for Fault Blocks 6-8. 
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TAB Aa ex 
WELL DATA—FAULT BLOCK 7 
Refer to Figure 14 


Well Rate,STB/Day Thickness Well Rate,STB/Day Thickness 





_No. (+Inj, -Prod) 1 oe No. (+Inj, -Prod) IDS. 
1 234 a9 13 -~282 47 
2 -264 44 14 ZSs 43 
5 234 39 US -288 48 
= -258 43 16 240 46 
5 234 bey, Nye -288 48 
6 -258 43 18 ZO 45 
7 226 38 9 SA 47 
8 -240 40 20 264 44 
2 Gee 2 21 -276 46 
10 -228 38 tag 258 43 
1] -270 45 Gs -270 45 
eZ Zac 42 24 258 43 


Production well rates during fill-up equal 25 STB/Day. 





TABLE XI 
WELL DATA—FAULT BLOCK 8 
Refer orem re. lt 


Well Rate,STB/Day Thickness 


No. (+Inj, -Prod) Ft. 
1 124 31 
2 -124 3] 
3 124 30 
4 -124 31 
5 -168 42 
6 168 42 
7 -168 42 
8 160 40 
y) -168 42 

10 156 39 


Production well oil rates during fill-up equal 15 STB/Day. 





Well 


INO; 


oN oO oH eR WwW NY 


Ne) 


10 
Wl 
12 
i> 


Production well oil rates during fill-up equal 15 STB/Day. 


Rate, STB/Day 
(+Inj, -Prod) 


-120 
-160 


-135 
-200 


-155 
-200 


TABLE XII 


Thickness 


ie: 
27.7 
30.4 
ag. 
Si 
30.4 
25.0 
53. 
40, 
3576 
40 
37.6 
40. 
40 


WELL DATA—FAULT BLOCK 9 


Peete LOn Hem ead 


Well 


No. 


14 
15 
16 
beg 
rs 
nS 
20 
ZA 
ae 
ae 
24 
ao 


Rate, STB/Day 
(+Inj, -Prod) 


-~150 


al, 
Pa 


-150 
-200 


-150 


-150 


-~150 


-150 


47 


Thickness 
iE. 


40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 
40 


* Program calculates injection well rates based on reservoir thickness 


and steady state material balance. 
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FAULT BLOCK 9 


IS ¥ | fee ke} 
1e ‘31 Jo? 03 
1a 6! Piogorns XA 
/9 /, 07 | 
{ oO 2 7 08 
| o — ‘dee od is / 
ee H 
x ° 2: hes ° 


Section 15 


FAULT BLOCK 10 


19 





Section 14 


xX Injection Well 
O Production Well 


Figure 15 Well Locations for Fault Blocks 9 and 10. 
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TABLE XIII 
WELL DATA—~—FAULT BLOCK 10 
INeLer toa toure 15 


Well Rate,STB/Day Thickness Well Rate,STB/Day Thickness 


_No. (*Inj, -Prod) Ft. No. (+Inj, -Prod) Ft. 
] -76 38 fo a 40 
2 s a6 22 -80 40 
5 ae 40 Zo -80 40 
+ -76 38 24 -80 40 
5 -70 55 Z5 * 40 
6 -80 40 26 -76 38 
7 -74 37 rag * 35 
8 35 Zo “ty 40 
9 a 40 29 -80 40 
10 -80 40 30 -78 ay, 
11 -74 oF eal -76 38 
il * 34 2 -72 36 
1 * 40 oD -70 35 
14 -150 40 34 *s eRe. 
15 -150 40 a5 -76 38 
16 * 40 36 a“ a7 
Ih -80 40 37 * 36 
18 -76 38 ae is ao 
19 -72 36 59 -66 oo 
20 i 34 


Production well oil rates during fill-up equal 8 STB/Day. 


* Program calculates injection well rates based on reservoir thickness 
and steady state material balance. 








APPENDIX B 


RATE ADJUSTMENT 
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ADJUSTMENTS OF INJECTION AND PRODUCTION RATES TO 
ACHIEVE UNIFORM WATER BREAKTHROUGH 
IN ALL FAULT BLOCKS 

When the model response was inspected section by section it was 
found that some of the sections were watering-out faster than others. 
This situation is also possible in the field. If the injection and production 
rates cannot be sufficiently adjusted, certain sections would be watered- 
out and shut-in while others would be run longer until ultimate recovery. 

If the rates can be adjusted enough the field can be operated so 
that all of the sections water-out at about the same time. The well rates 
used in the model were adjusted, with the total field injection rate held 
constant, to determine the effect of uniform water breakthrough from 
section to section. 

The results, as shown in Figure 16, indicate that for the same 
period of operation the oil production was increased and the water 


production was decreased. 
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APPENDIX C 


COMPUTER PROGRAM 
iW) FP rogramelbisting 
2, Sample Data Deck 





(DAO N|N 20:90 N80 A'A0 A000: A0N0ABAAAN1ANDN AMA AA ANA AD AAA AAAAAAAAAADAANDAAAGGAAAAD 


PROGRAM OL8NKVT( INPUT, OUTPUT) 
A STREAMLINE MODEL FOR WATERFLOOD SIMULATION 


THIS MODEL IS SPECIFICALLY DESIGNED To SIMULATE 4A 
WATERFLOOD OPERATION IN WHICH AN OIL BANK IS FORMED AND 
MOVES THROUGH THE: RESERVOIRe THE MODEL 15 CAPABLE OF 
HANDLING RESERVOIRS WITH THICKNESS VARIATIONS ANO: 
IRREGULAR BOUNDARIES. 


REFERENCES? 


LEBLANCsJelLes AND CAUDLE®B.He, A STREAMLINE MODEL FOR 
SECONDARY RECOVERYs SOC, PET. ENGRe Jo (1971) II 97~ 
RUSTeCeoBes A STREAMLINE MODEL FOR OIL BANK BUILEDIP IN 
A WATERFLOODs M.Se THESIS, THE UNIVERSITY OF TEXAS AT 
AUSTIN (MAYs 1972). 
WESSELSoeJeWes APPLICATION OF THE STREAMLINE RESERVOIR 
MODEL TO RESERVOIRS HAVING VARIATIONS IN FORMATION 
THICKNESS® MeS. THESISs THE UNIVERSITY OF TEXAS AT AUSTIN 


DOCUMENTED BYt KeE.GOLTZ (MARCH,1975) 


SFKSSOROHGHaGonen eat huatdavyaeteatoua 


® SCHEMATIC PROGRAM FLOW CHART & 
SHeHesovtastonnartettovtdavpaveotouag 


I READ INPUT DATA I 


I 
I COMPUTE STREAMLINE STARTINS POINTS, I 
I SOURCE AND SINK PRESSURES I 


I 


cteetenoerecnenworowey §6JTERATE EACH TIME PERIOO I 


O@ 0-6 P=4 0-6 Og D-4 0-6 D=4 O-g 4 Og D-4 Oe OG OG 0-4 Dt OO 


I 
wrcwwetowerororo!] §=6JTERATE EACH PRODUCTION WELL I 
I Se®@ese eee GCoeee@goeZFteaoaes@Fao ®®a ea g@tF®oFe ¢o 
I I 
I SP@eeeeeSE@Pgoe®aauae Se@tPFaas®a Gyeo*® a GG 
T eworecccocwcorrey] §=ITERATE EACH STREAMLINE I 
I I Seee Sees Gvuaeeaave SEP ®an Fe © geteug Fo 
i 1 I 
I I SSSSF OCS COeseCeFeSOge®Gaweg SSS SF GeO get FF eFoaaaeGda ®& 
I JI eecery COMPUTE VELOCITIES AND TIME INCREMENTS I 
ie) ie: I AT EACH POINT ALONG STREAMLINE I 
J J J Se eF®@ue S®eeeeeesPGgeFeucas VZePFFaeGe@ geoFGuat*eoauatea® 
p ie I 
J I J Pees ee SeCeeesSSCeeFqgoF®eaes® ©S8SFeeFe Gg eZFaesonvea SF G®ee' 
I 1 I I MOVE FLUIO PARTICLE ALONG STREAMLINE USING I 
I 1 it I FINITE DIFFERENCE APPROX. OF DARCYS LAW I 
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inary I ACCUMULATE STREAMLINE PRODUCTION 1 

I I Dk ek eee eee ee St oP eee Ce ee ST 

ieee ae I 

I I I OOP CES SO CSCS GB e® Ge we FTF ®% HOF e BESTEST CAG. 

| AS) Gee I TEST FOR BREAKTHROUGH OF FLUID I 

I I iT I PARTICLES INTO PRODUCTION WELLSe I 

ae I ADJUST STREAMLINE PRODUCTION t 

I I I I RATES AS REQUIRED I 

I I I 088 88 8888 68S ww FO wee SOSVF% GSES © ZETVSE* DOOD 

I I I I 

I I I NO S®e e@ SF SEF geo®awoen Ze®® gaFa &get*®Gn Fe ee 

I IT seecewrrenwnw] TEST FOR END oF TIME PERIOD I 
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I I NO O88 S22 0% we FW wUe se S%enaFa GyeF2u®aa- 

T ewerweeeewenrel «TEST FOR ALL STREAMLINES I 

I Se ee @@2eaeae@ge “ava awe PFoaPVP® ea Fa Gy st Sea 2a: 

I I YES 

I NO ede ke kK de ed ee ee Le ee ee ee ee 

wer weet ouawuncee] TEST FOR ALL PRODUCTION WELLS I 
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I PRINT/PLOT WELL ANO FIELD RESULTS I 
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CV SWFVF ww BV aC eZ eZ eee *u~u | TEST FOR TOTAL TIME I 
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INPUT DATA DEFINITIONS 


POR = POROSITY (DECIMAL) 
SCw = SATURATION CONNATE WATER (DECIMAL) 
SOI = SATURATION OIL, INITIAL (OECYMAL) 
SOR = SATURATION GIL RESIOVAL (DECIMAL) 
SGR = SATURATION GAS RESIOUAL (DECIMAL) 
SGOR = SOLUTION GAS#OIL RATIO (SCF/BBL AT RESe CONDITIONS) 
NP = NUMBER OF PRODUCTION WELLS 
NI = NUMBER OF INJECTION WELLS 
NBW © NUMBER OF IMAGE WELLS USED TO BOUND PATTERN 
NBP = NUMBER OF POINTS USED TO DEFINE BOUNDARY 
NST = NUMBER OF STREAMLINES FOR LOWEST RATE PRODUCTION WELL 
NXH © NUMBER OF X@DIRECTION ENTRIES IN THICKNESS GRID 
NYH = NUMBER OF YeDIRECTION ENTRIES IN THICKNESS GRID 
NEOIT = 0 OR 1 
IF 9 OR BLANK@PROGRAM WILL BALANCE PROD/INJ RATES 
IF 3 PROGRAM INTERNAL BALANCE FEATURE IS DEFEATEO 
H ® AVERAGE RESRVOIR THICKNESSsFT (USED FOR IMAGE WELLS) 
RI = RADIUS OF STARTING AND ENDING CIRCLESo FT 
DELT =» TIME INCREMENT (DAYS) FOR UPDATING CONSTANT PRESSURE 
CORRECTION FACTOR» SINCE DELT MUST START OUT SMALL BUT 
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62 
63 
$4. 


65: 
67 
639. 
69: 
70 
7} 
72. 
73 
74 
75 


77 
78 
793 
89 
3} 
82 
B3 
B4 
85' 
95. 
87 
88 
B63 
990 
9} 
92 
93 
94 
95' 
95: 
97 
99. 
99 
100 
101 
102 
103 
104. 
lo5' 
105: 
107 
199 
109 
lio 
lil 
li2. 
113 
lis. 
115: 
115: 
$17 
113 
11> 
120 
21 
122. 
123 





NONAANANAVANAGAA AAAADAXANDAAAQ 


aoaaqn 


BECOMES LESS CRITICAL AS THE SOLUTION PROGRESSES, 

THE VALUE OF! DELT IS MULTIPIED BY TEN AFTER THE 

OIL BANK BREAKS THROUGH IN THE FIRST STREAMLINES 
DELP = TIME INCREMENT FOR PRINT AND PLOT (JAYS) 

FIRST PRINT/PLOT IS MADE AT DELTsTAHEN EVERY DELP: 
TMX = MAXIMUM TIME FOR SOLUTION (DAYS) 
QMNP = LOWEST PROOJCTION WELL RaTe (STB/DAY) 
QMNI = LOWEST INJECTION WELL RATE (STB/DAY/FT) 
XMX ws ABSOLUTE X OIMENSION AT WHICH STREAMLINES AR= CUT OFF: 
YMX s& ABSOLUTE Y DIMENSION AT WHICH STREAMLINES ARE] CUT OFF) 
GMBsOMSewMB = GAS*eOIL AND WATER MOBILITIES 
XW(I),YWC(I) = XeY COORDINATES OF wELL LOCATION 
HI(I) = THICKNESS OF RESERVOIR AT WELLoFT 
Q(I) = WELL INITIAL: OIL PRODUCTION RATE (ST3/DAY) 
QO(I) = WELL PRESENT (FILL=UP) OIL PRODUCTION RATE (STB/0AY) 
GOR(I) = WELL GAS@OIL RATIO (SCF/SB8L AT RES. CONDITIONS) 
XB(I)sYB(I) = Xe¥ COORDINATES OF BOUNDARY POINTS 
XH(I),YH(T) = Xe¥Y COORDINATE VALUES FOR THICKNESS GRID 
HT(IeJ) = VALUES OF! THICKNESS CORRESPONDING TO COORDINATES 
OF THICKNESS GRIO 


FHPOSHGHEEHHHOHOHHHHEOHHHOGHOTHEGO KE EPHHEEESHOHEEHHe Heese eeseegeshieagad & 


COMMON/AAA/S xW(125) eo YW(125) op HI (125) 01125) 0X8 (125) » 
1Y8(125) »sNBPeNPoNI eo NBWoet 

COMMON/BBB/ €C(55955) 2D(55) 9Q1 (55) 

COMMON/CCC/ X(80020) s¥ (80920) 9X0 (80920) 2¥u(80920) » 
1WQ(80».20) oSCL oXMXeYMXeNST se Q@MNPsRGO (80920) » 
2SGOReNPIWe TMETI 

DIMENSION TSUM(80020) PP (80520) oP] (80920) 209 (80920) » 
1CW(80 420) »Q0(125) 950R(125) 909(80220) » 
2XH(20) oe YH(20) osHT (20920) 

READ 530.5 PORsSCW95012SOReSGR#SGOR 

READ S400 NPoNIoNSagNBPeNST 

READ 540> NXHeNYHe NEDIT 

READ 53059 HeRIsOELTsDELPs TMX» QMNP 

READ S309 QMNIoXMXeYMXeGMB2 ONS eWMB 

DELTT=DELT*10. 

NCOUNT#0 

PLTFT=009*®XMX/YMX*0 445 

TOTFTs0e7¢PLTFT 

NPI WeNPen] 

NPIBw=NPIWeNBw 

READ S509 (XW(I) oYA(T) pHI (1) 5 Q(T) ,Q0(T) pGOR(T) sI=1,NPIBW) 

READ S609 (XB(1) eY3(1) sIT=19N3P) 


CONVERT SURFACE RATES TO RESERVOIR RATES 
FVFle le2 


D0 9 Jal ,NPIW 

Q(I) 2Q(1)#1le2 

Q0(1)2#90(1)#1.2 

QMNPEQUNP®@le2 

QMNIEQUNI #102 

PRINT 590 

PRINT 6808 NPeNIeN3weNBPoNST 
PRINT 6900 HeDELT oe ITMX,RI 

PRINT 5002 DELPe®XMXeYMX»QmMNP 
PRINT 5100 POReSCweSOleSOR®SSR 
PRINT 520° OMBeWMBeGMB 

PRINT S500 (XWCI) eo ¥W(I) sHI (I) 0 Q(T) 2QO0(1) eGOR(T) eI 81 eNPIBW) 
PRINT 5300 (XB(1) eY8(I) slelsNBP) 
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1246. 
123' 
B26, 
127 
128. 
1293 
130 
13) 
132 
133 
134. 
133: 
135. 
137 
139 
139 
149 
1461 
142. 
143 
166. 
145: 
1465, 
147 
143 
149 
150 
15) 


is 3’ 


17} 
172. 
173 
176 
173! 
17s 
177 
173 
1793 
185 
18) 
182. 
383° 
184. 
185: 





Oa00 


10 
20 


30 


31 
3¢ 


(eR ek @) 


40 


aaqg 


50 


60 


70 


80 


90 


200 


BALANCE INJECTION WELL RATES CONSIDERING RESERVDIR 
THICKNESS ANDO STEADY STATE OPERATION 


IF(NEDIT oGT. 0) GO TO 32 
SMQI29.0 
SMQP329 20 
DO 20 I=l]esNPIW 

IF (Q(1)+LTe0.0) GO TO lo 
SM@IeSMQT+Q(7J) 

Go 10 20 
CONTINUE 
SMQPESYQP4+Q(T)/HI (I) 
CONTINUE 
00 30 I=leNPIW 

IF (Q(T) e¢L T2020) 60 TO 30 
Q(1) 209(1) @SMQP#HI (1) /SMQI 
CONTINUE 
Go TO 32 
PRINT 595 
CONTINUE 


CHECK IF OIL BANK WILL FORM 


IF (SOReLTeSOI) GO TO 40 
PRINT 570 
GO 10 +70 
CONTINJE 


CONVERT FIELD UNITS TO CGS 


00 50 J=leNPIBW 
XW(T) aXW(1T) #30248 
YW(T)ayW(1)%39.48 
Q0 (1) =90(1) #1.84%01 
HI (Tt) eH1(1)*30.48 
Q(T) 2Q(1) 41.840) 

IF (I,LEeNPIW) SO TOD 59 
HI (1) 2H*30.48 
CONTINVE 
DO 60 IT=]¢eNBP 
XB (7) aXB(1) #30248 
YB(I) syB(1)#30,48 
CONTINUE 
Heh*30.48 
DO 70 I=leNXH 
READ S805 (HT (IeJ) e Je] eNYH) 
CONTINVE 
READ SS90¢0 (XH(I) oe TS1eNXH) 
READ 580+ (YH(I) ¢I31eNYH) 
DO 89 141 eNXH 
00 89 J=leNYH 
HT (Le J) SHT (19 J) 930048 
CONTINUE 
00 90 [8] ¢eNXH 
XH(I) e&XH(I)#390.48 
CONTINUE 
DO loo T2leNYH 
YH(1) aYH(T) #30.48 
CONTINUE 
RIERTS30.48 
DELTEDELT#86400- 
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136, 
187 
199. 
183: 
190 
191 
192 
193 
134. 
195! 
195, 
197 
139 
139 
C00 
201 
e02. 
e03 
204. 
e05' 


e07 


e093 
e10 
eli 
e12. 


e25' 


223 
223 
230 
e3} 
¢32 
e33 
234. 
235' 
235: 
e37 
239. 
e393 
240 
24) 
242. 
2463 
C44, 
245) 
CH: 
247 
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QNooO|0 20 


OOOO 


110 


DO.0O0 ODAADAANOO0 A8AAAAAHAADADAIA AAD AAAGDAAAGAGAAD 


OELTTmDELTT*86400. 
DELPSDELP#86400. 
TMX=TMX®86400¢ 
QMNP2Q4NP#].8401 
XMKBXMX#30048 
YMX=YMX®30.48 
QMNI=QUNI#1.8401 


SUSROUTINE BDUND CALCULATES THE NECESARY RATES 
FOR PREVIOUSLY SELECTED IMAGE WELLS SO THAT THE 
STREAMLINES ARE CONFINED WITHIN THE RESERVOIR 
BOUNDARIES, 


CALL BOUND 
PRINT 600 


CHECK QMNI AND CALCULATE VELOCITIES NEAR 
INJECTION AND PROOUCTION WELLS 


DO 119 I=l»NPIBW 

IF (Q(1)eL T2020) GO TO 110 
IF (Q(I)/HI(I1) eL Te OMNI) QMNIEQ(T) sHI (I) 
CONTINUE 
VLCT=9.5°QMNI/RI 
VLCP=0,4#QMNP/H/RI 
NTSLeQ 


THIS SECTION STARTS A STREAMLINE ON THE CIRCLE OF A 
PRODUCTION WELL AND FOLLOWS IT BACKWARDS TO LOCATE THE 
PROPER INJECTION WELL. THIS PROCESS IS CARRIED OUT FOR 
ALL PRODUCTION WELLS AND ALL STREAMLINES» THE 
STREAMLINE STOPS AT A POINT WITHIN THE CIRCLE OF AN 
INJECTION WELL. THIS POINT BECOMES THE STARTING POINT 
WHEN THE STREAMLINE IS FOLLOWED IN THE FORWARD OIRECTION 
TO ACCUMULATE THE: PRODUCTION HISTORY. 

AS THE STREAMLINES ARE BEING FOLLOWED A PLOT OF) THEIR 
TRACE IS MADE>s (NOTE THAT ALL PLOT STATEMENTS HAVE 
ASTERISKS In THEIR STATEMENT NUMBERS») 


SSHSSEHHEEHEHEHEHHEHHFEHEFECEO HE HOHSZOHEHEH HERES HATH EHEHHOOEEC HHH & 


PROGRAM VARIABLES 


J 8 INDEX FOR WELLS 
K ® INDEX FOR STREAMLINES 


XIsV¥I = POSITION ON STREAMLINE OF FLUID PARTICLE BEING TRACED 

VXeVY = VELOCITY OF: FLUID PARTICLE IN X AND Y OYRECTIONS 

PP © PRESSURE AT PRODUCTION WELL (SINK) 

PI] = PRESSURE AT INJECTION WELL (SOURCE) 

X(JoK) eV (JeK) = STARTING POINTS FOR WATER FRONT FRI INJe WELL 
XO(JoK) eVYO(UeK) = STARTING POINTS FOR OIL FRONT FROM INJe WELL 
TSUM(UeK) = SUMMATION OF TIME 

RGO(JU,«K) & STREAMLINE GASeOIL RATIO 

OQ(JeK) = STREAMLINE OIL RATE 

CO(JeK) = STREAMLINE CUMMULATIVE OIL PRODUCTION 

WO(JeK) & STREAMLINE WATER RATE 

CwiJeK) & STREAMLINE CUMMULATIVE el PRODJCTION 

SCL # PLOT SCALE FACTOR 

NSL # NUMBER OF STREAMLINES FOR THE PRODUCTION WELL 


58 


243. 


e55' 
255: 


264. 
2635! 
C56: 
267 
263 
269: 
2709 
27} 
272 
273 
274. 
2735! 
275. 
C77 
279. 
e793 
230 
28) 
282, 
233 
#ategy. 
#A#2832) 
#ate2gs,. 
287 
289. 
2893: 
290 
291 
292. 
293 
294. 
295! 
295. 
297 
259, 
299° 
309 
301 
302 
303 
3046. 
305' 
306: 
307 
#A°3939. 

A 309 
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OOCORVRN0 


Onoan 


129 CP CF O'OD OC 


120 


130 
140 


150 
160 


170 


BEING WORKED ON 
OT = TIME INCREMENT FOR SMALL STEP TAKEN. STEP SIZE IS 
EQUAL TO e254R], RI NEEO NOT BE SMALLER THAN 
1/50TH OF THE DISTANCE BETWEEN WELLS. 


SHSeshGe est SGHOGGGsGHGhaessessets taste gee SGGGnagsedtpHeeGnGGenansedeGnann ag 


SCLS5,238/YMX 
XME2.257XMXPSCL 
CALL BSNPLT (4LPLOTsXM925925) 
CALL PLT (1el25*#XMX#SCLo1.125*%YMX#SCL9=3) 
CALL PLT (XB(NBP) *SCLsYB(NBP) *SCL,3) 
D0 129 IT=l»sN3P 
CALL PLT (X8(I)*SCLe YB(1) #SCL92) 
CONTINUE 
DO 149 ImlsNPIW 
IF (Q(1)eGE.0e0) GO TO 130 
CALL SYMBOL (XW(I) *SCLoYW(I) #SCL9,07919090)1) 
GO TO 140 
CONTINUE 
CALL SYMBOL 
CONTINUE 
DO 229 J=leNPIw 
IF (Q(J)eGE20.0) GO TO 220 
NSLEeNST#Q (J) /QMNP 40,5 


(XWCT)*®SCLeYW(1) *SCLs,0797009)) 


THE STREAMLINE STARTING POSITIONS ARE EVENLY SPACED 
AROUND THE PRODOCTION WELL 


DO 210 Kz=lsNSL 
MI eXw( J) oRI@COS (2 20°K#6,2832/NSL) 
YIZYwW (J) eRI@®SIN(1 e00K86,2832/NSL) 
CALL PLT (XT #SCLsYI#SCL 33) 
PP(J9K) =020 
DO 159 L=l»sNPIS8w 
PP(J9K) =PP(JeK) SQ(L) SHI (L) *ALOG( (xTeXW(L)) FH2¢(VYIeVw(L)) ##2) 
CONTINVJE 
CONTINUE 
IF (ABS (XI) -GTeXMX) 
IF (ABS(YI) -<GTeYMX) 
VX20.0 
VY=0.0 
DO 179 L2alsNPI BW 
DEN=(XIT=Xw(L)) ##2*(YIoYwi(L) ) #2 
VX=EVXeQ(L) (XT eXwW( lL) )/(DEN@HI (L) ) 
VYEVVeQ(LIF(VIOVYW (LL) S(DENPHI (L) ) 
CONTINUE 
VTESQRT (VXF*#oeVYee2) 
OT20.25*RI SVT 
XTEXTe+DT@vx 
VYIRSY]e3T avy 
CALL PLT (XI#SCLeYI*SCL 92) 


GO To 180 
GO To 180 


CHECK IF THE STREAMLINE HAS REACHED AN INJECTION WELLe 
THE VELOCITY OF THE FLUIO PARTICLE IS CHECKED FIRST TO 
SEE IF IT IS LARSE ENOUGH TO BE IN THE VICINITY OF AN 
INJECTION WELL. IF IT IS_y EACH INJECTION WELL IS THEN CHECKED. 


IF (VTeLE.VLCI) GO TO 160 
DO 200 LL&lsNPIBw 
IF (Q(LL).LE«000) GO TO 200 


2), 


310 
3}] 
312, 
313 
316. 
313) 
315: 
317 
@A%3193 
#A€#379 
PA%I2H 
@A%32) 
#A8322 
#A%323 
#Ae324 
@A#325) 
@#A*325, 
#4%327 
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A323 


#A%329 
#A%330 
#a%33) 
#A®332. 
333 
334 
335: 
335 
337 
333 
339 
349 
34} 
34? 
#343 
34%. 
345! 
34%. 
347 
343 
349 
350 
35} 
352 
353 
354. 
3535’ 
355. 
357 
359 
359 
380 
36} 
#362 
363 
364 
365 
365. 
367 
363. 
369 
370 
371 
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180 


190 


200 
210 
220 


230 


RADSSQRT ( (XT eXWILL) ) ##2¢ (YI mYW (LL) ) #2) 


IF (RAD.GTeRI) SO TO 200 
CONTINUE 
NTSLENTSL4] 
X(J9K) XI 
VY(UeK) xYI 
TSUM(JeK) £020 
XO(JeK) EXI 
YO (JeK) SYI 
PI (JeK) 2020 
DO 199 LK=leNPIBW 


PI (JeK) SPI (eK) @O( LK) SHI (LK) MAL OG ( (XT @XWILK) ) #24 (YT eYW(LK) ) @@2) 


CONTINUE 
RGO(J,K) =GOR (J) 


IF (SCW*SOI*SGReGEel.0) RGO( JK) 25GOR 


OQ(JekK) =Q0(J)/NSL 
CO(JeK) =0-0 
WO(JeK) 20290 

GO TO 210 
CONTINUE 

GO TO 160 
CONTINUE 
CWlJeK) =020 
CONTINUE 
CALL PLT (00020009999) 
PRINT 610 


END OF RUNBACK SECTION 


CTR®1.,0 
PRINT 590 
TSTOP=HDELT 
ISMP=DELT 
NBPl=N3P¢} 
CONTINUE 
SMCTYx000 


THIS SECTION STARTS THE STREAMLINE AT THE INJECTION 
WELL FOUND ABOVE ANDO FOLLOWS IT IN THE FORWARD OJRECTION 


TOWARD THE PRODUCTION WELL, 
WATER FRONTS ARE CLOSELY wATCHED 
BREAKTHROUGH INTO THE PRODUCTION 
HISTORIES ARE ACCUMULATED AS THE 

LOOPS 420 AND 410 ARE DONE FOR 
STREAMLINES UNTIL: TIME EQUALS TSTOPe. 
PRESSURE CORRECTION FACTOR 
RESULTS ARE PRINTED IF DESIRED. 


MOVEMENTS OF TH OIL AND 
TO DETERMINE WHEN THEY 
WELL. THE PRODUCTION 
FRONTS MOVE. 

ALL WELLS AND ALL 

THEN THE CONSTANT 
IS UPDATED AND THE 


Seta GetegaeteGaetaoaeseGedhaggaetcectaatogaga dette ceagn,eSGaotsgrangdeGetany 


ADOITIONAL: PROGRAM VARIABLES 


CTY = CONDUCTIVITY RATIO 


CTR = CONDUCTIVITY RATIO CORRECTION FACTOR 


TSTOP = TIME TO STOP AND UPDATE CTR 
TSMP &@ TIME TO PRINT ANO PLOT 


XIeYI =# LOCATION OF WATER FRONT ON STREAMLINE BEING FOLLOWED 
XOlsYOI = LOCATION OF OIL FRONT ON STREAMLINE BEING FOLLOWED 


FESSHEHHS SHEESH EEEEESHEFEKETOSETHRGEEHOTHAEEEHAEHHHHHEHETHTAOG HY 
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60 


372. 
373 
376. 
375' 
375. 
377 
373 
373 
389 
38) 
382. 
383 
384. 
385' 
385 
387 
389. 
383 
399 
39] 
392 
393 
39%. 


#395, 
3397 
393 
399 
600 
601 


§93 
$04 
695) 
605% 
407 
$93 
$093 
$19 
$1) 
&}]> 
413 
414. 
615° 
415, 
617 
$19. 
413 
$20 
$2) 
422 
$23 
$24. 
G25! 
626 
627 
423 
$2) 
630 
$3) 
630 
433 





Teli ele) 


oOoQdan 


240 


250 
260 


270 
280 


290 
300 


310 
320 


330 


429 JaloNPIW 
a Ten On Ds GO TO 420 
NSLEeNST#Q (JU) QMNP 90,5 
DO 4190 KeloeNSL 
zs 
PAA ee tank) «GT. Kus GO To 410 
IF (ABS(Y(JeK))eGT,YMX) GO TO 410 
IF (WO(JU9K) eGTe000) GO TO 409 
XJTBX( JK) 
YTSV¥ (eK) 
TSUMI sTSUM(JeK) 
XOTEXO(J9K) 
YOIT=YO(JoK) 
CONTINUE 


DETERMINE THICKNESS OF RESERVOIR AT LOCATION OF FRONTS 


00 250 M=leNXH 
IF (XTeLT.XH(M)) GO TO 260 
CONTINUE 
MaNXw 
CONTINVE 
DO 279 NzleNYH 
IF (YIeLTeYH(N)) GO TO 280 
CONTINUE 
NaNYH 
CONTINUE 
HTP=HT (MeN) 
Vee (RGO (JyK) eLE*SGOR) GO TO 369 
DO 299 MzleNXh 
IF (XOTeLTeXH(M)) GO TO 300 
CONTINUE 
M=NXH 
CONTINYE 
DO 319 NeleNYH 
IF (YOIeLTeYH(N)) GO TO 320 
CONTINYE 
N2NYH 
CONTINUE 
HYPOsHT (MeN) 


TS AND. ACCUMULATES 
ON MOVES THE OIL AND WATER FRON 
PRODUCTION BEFORE. THE OIL FRONT REACHES & PRODUCTION WELL, 


VX0=0.0 
VY0=0.0 
PO20.0 
VX20.0 
VY=0.0 
eee: l»NPIBW 
=1°9 
S00= (KOT=AW(L)) #820 (YoL“¥w(L}) #42 
EPO@9(L) /HI(L) *ALOG(SQ 
Oe ral ts bar or-Xw tL.) 17800 
VYORVYO*Q(L) HI (L) ®(YOT=YW(L))/SQ0 
SQE(XIToAw(L) ) P20 (VTeYW(L) ) #ee2 
PePeQ(L)/HI(L) #AL03(SQ) 
VX3VX09(L)/HI(L) (XT eXW(L)) 7SQ 
VYRVYOO(L) SHI (L) @ (YI eYw(l)) 7$Q 
oe eat K)) PO=xPP(JeK) J 
7 = 
GHGs) (cot 10k) “PONG + (ROEP) 0M36(P PI(JeK)) 
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61 


634. 
$35' 
635. 
$37 
$33. 
$39 
649 


645) 


46} 


$35! 


667 
668 
$69 
$70 
47) 
472. 
473 
474 
4735) 
975: 
477 
473 
473 
$89 
$3) 
682. 
$93 
684. 


685s, 
687 
$993 
49) 
$99 
$9) 
692 
693 
694. 
$93! 





aaanoaaon 


349 


‘Oao oO 


346 
350 


ODog 


360 


370 


WMB) ®GMB) 
VLXEVX/ (60 2B3Z2*POR® (LeSCWHESOR=SGR) D#CTY/CTREHI (UI SATP 
VLYEVY/(6e28328POR*( 1LeSCWHeSOR@SGR) )ACTY/CTROHMI (J) STP 
VELESORT (VLX@e2eVLY Heo) 
OTBRI/VEL 
CRFE1,0*0Q( Je K) ®NSLECTRS(Q(U) *CTY) 
VLXOEVXO0/ (6. 2832#POR#( 1 eSCWeSOTeSGR)) HCTY/CTRECRE SHI (J) /HTPO 
VLYO=VYO/ (602832#PIR® (1 -SCWeSOT@SGR) ) FCTY/CTRECRE Hy (J) /NTPO 
VELO2SQRT(VLXO08*2¢VLYO#%2) 
IF (VELO.GTeVEL) OT=RI/VELO 
IF (TSUMI*¢OT.GTeTSTOP) OTaTSTOPeTSUMI 
XTEXToOT#VLX 
YISVIOOT#VLY 
TSUMIZ=TSUMI*DOT 
XOT=XOI*OT#VLXO 
YOTSYOI*OT#VLYO 


CHECK IF OIL FRONT HAS REACHED A PRODUCTION WELL. 
THE WATER FRONT IS ALSO CHECKED JUST IN CASE IT 44S GOTTEN 
AHEAD OF THE OIL FRONTe THIS OF COURSE IS NOT NORMAL 
BUT CAN HAPPEN IFi THE OIL PRODUCTION RATE DURING 
FILLeUP IS TO LARGE. 


IF (SQRT((XIT@Xw(J)) 20 (VYToYw(g)) ##2) eLEeRI) GO TO 390 
CO(JeK) =CO(S9K) *OT#OQ( eK) 
IF(VELD elEe VLCP) GO TO 350 
IF (SQRT((XOLeXW(S)) 09206 (YOT@¥W(J)) #2) eLEeRI) GO T): 348 
DO 369 LX=leNPIBw 
IF(Q(LX) eGEe 000) GO TO 349 
IF (SQRT((XOTeXWILX) ) @@26(YOTYW(LX)) #2) eLEeRI) GO TO 348 
CONTINJE 


WHEN THE OIL FRONT REACHES A PROOUCTIOV WELL THE 
STREAMLINE GAS*OIL RATIO IS SET EQUAL TO SGOR, 


GO TO 350 
RGO(J,K) =SGOR 
NCOUNTS] © 
CONTINVE 
IF (TSUMI,LTeTSTOP) GO To 249 
X(J9K) 2X] 
Y(JeK) e&YI 
XO(J9Kk) =XOI] 
YO( JK) =YO!I 
TSUM( JK) ETSUMI 
GO TO 620 


THIS SECTION MOVES THE WATER FRONT AND aCCUMULATES 
THE PRODUCTION AFTER OIL FRONT QGREAKTHROUSH 


CONTINUE 

VX®0.0 

VY=0.0 

Pe0e0 

DO 379 Lel»NPIaw 

SQB(XyeXwihL) )##2e(VJoeyw(lL)) ase 
PzPeQ(L)/HI(L) #ALOS(SQ) 
VXEVXOO(L)SHI(L) @(KyeXW(L)) 7S 
VVEVYeQ(LISHICLI A(T eYW(lL))/S@ 
CONTINUJE 

IF (PeLTePPl(JeK)) PePP(J9K) 
CTYS(PP( eK) @oPI(J9K))/(( (PP (9K) op) SOMBe (Pedy (JeKk) ) /WMB) ®GMB) 
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333 
334. 
335: 


337 


564. 


3463 
549 
559 
351 


353 
254. 
355: 
355: 
357 
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380 


390 


'CHOVOO 


400 


ocd 


410 
420 


Ooaqaaadaa 


VLXSVX/ (6, 28328POR*( LeSCWHSOR=SGR) ) OCTY/CTRENI (UU) SATP 
VLYEVY/ (6e2B32#P0R® (LeSCWHSOR@SGR) )MCTY/CTROEMI (UU) Z4TP 
VELESQRT (VL XeS2¢VLYee@2) . 

OTSRI/VEL 

IF (TSUMI¢DT.GTeTSTOP) DT#=TSTOP=TSUMI 

XIEXTeDT#eVLX 

YISYIeOTeVLY 

TSUMI&TSUMI*OT 

OQ (Je K) ==Q(J) OCTY/ (CTRENSL) 
CO(U9K) =CO( Ue K) HOT #0Q( Je K) 


CHECK FOR WATER FRONT BREAKTHROUGH 


IF (SQRT((XIT@XW(S)) #820 (YT e¥Ww(y)) #2) eLE.RI) GO TD 390 
CONTINVE 
IF (TSUMI,LTeTSTOP) GO To 249 
X(J9K) eX] 
Y(JoK) YI 
TSUM(JeK) =TSUMI 
XO(UeK) =XOI 
YO(JekK) =YO!I 
Go TO 410 
CONTINUE 


THIS SECTION IS USED WHEN THE WATER FRONT FIRST 
BREAKS THROUGH, 


RGO(JUeK) =SGOR 

CTY=wM8/G6MB 

WQ(JeK) =e Q(J) OCTY/S(NSL®CTR) 

OQ (JK) =020 

X(JoK) 2X] 

Y(JeK) YI 

TSUM(JeK) =TSUMI 

XO (J—9K) =XOI 

YO(J9K) =YOI 

CWl eK) B(TSTOP@TSUMI) OWO(JeK) 
GO TO 410 

CONTINUE 


ACCUMULATE PRODUCTION AFTER WATER BREAKTHROUGH 


CTY=wM3/GMB 
WO(UeK) =9Q(J) OCTY/ (NSL*®CTR) 2 
CW (Sok) HCW (Jo K) *DELT OW (Uy K) 
SMCTY2SMCTY*CTY 
CONTINUE 
CTRESMCTY/NTSL 

IF (TSTOPeLTeTSYP) GO TO 450 


AT SPECIFIED TIME INTERVALS 

ALL STREAMLINES ARE SUMMED AND RESULTS ARE PRINTED 

FOR EACH WELL AND: FOR THE ENTIRE FIELDe 

PLOTS ARE MADE OFI THE OIL FRONT AND WATER FRONT POSITIONS, 


ORT20,0 
WRT=0.0 
SOIL29.0 
SWATs9.0 
DO 449 JaleNPIw 
IF (Q(J)}eGE.20.0) GO TO 440 
GRE0.0 


@ 
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559. 
353 
560 
2$) 
562 
363 
566. 
565 
365. 
367 
363 
3693 
270 
371 
372. 
373 
3746 
375' 


593) 


2935 
399 
600 
60) 
602 
203 
6804. 
605) 
60S: 
607 
503 
$093 
210 
#511 
512 
$13 
514. 
615) 
615. 
Se 
6395 
613 





OOo ago 


aod annann 


430 


440 


441 


450 


460 


470 
480 


490 
500 
$10 
$20 


OR®0.0 

WR80.0 

OPRDEO,.0 

WPRD209,0 
NSLE@NST#Q (JU) /QMNP*0,S 

00 439 KeleNSL 
OPRD=OPRD*CO(JsK)/158990- 
WPRD=WPROD*CW(JeK) 41589906 
GREGReRGO(U9K) ©0Q(JeK)/1 28491 
ORSO0R40Q(J9K)/1 08401 
WREWROWQ(J9K)/108401 
CONTINUE 

GSR=GR/0R 

WOR=WR/OR 

EX#Xw(J)/30.468 
WYEYW(J)/30048 

PRINT 6209 JoEXowy 

PRINT 6309 GReORpWR: 
PRINT 6409 OPRDewWPRDsGSReWOR 
ORT=ORT*OR 

WRT=wRT¢wWR 

SOILsSOIL*¢OPRD 
SWATESSWAT¢AWPRD 

CONTINUE 
TMEI=TSTOP/86400.0 

PRINT 65059 ORT sWRT 

PRINT 651s SOIL sSWAT 
PRINT 6609 TMEIeCTR 


MAXIMUM PLOT LENGTH IS 11 FT. 
PLOT LENGTH IS ESTIMATED ANDO IF IT IS GREATER THAN 
11 FT» THE PLOTTING IS STOPPED. 


IF((11,0"TOTFT) eL Te PLTFT) GO TO 44) 
TOTFTSTOTFT¢PLTFT 

CALL O3PLT 

CONTINJE 

TSMP=TSMP*DELP 

CONTINUE 


REMENT TIME TO NEXT STOPPING POINT. 
fee IME IS GREATER THAN TMX = STOP CALCULATIONSs | an 
IF OIL FRONT HAS BROKEN THROUGH IN AT LEAST ONE STREAMLINE 
sNCOUNTELSSET OELT = DELT#10, 


IF (NCOUNT.EGe1) DELTSOELTT 
TSTOPaTSTOP*DELT 
IF (TSTOP.GTeTMX) GO TO 460 
Go TO 230 
CONTINUE 
CALL ENOPLT 
Sa wiicek rene Se el39* NI c S9I3e* NIM = SeT30% NBP cg 913 
®o13) : 
ee ye ZS 9F10.59% DELT 2 *eFLO.S5e* THX = oF loels 
EARN 7? 48% + *OELP © eF1 90590 XMX © HeF1 Hel s*® YX 2 FoF Ogle 
@#,F10,5) 
a) Bee oF 6eSe® SCW & oF 659% SOL = Meh beS5ee SO 
#,F6 #® SGR s @F6.5) 
aor e ee soi woe 2 #,F10,5,%: WAT MOB = #9F 10059 GAS MOF 
ls #9F 1005) 
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220 
62) 
$22. 
623 
524%. 
$23! 
626. 
627 
623 
929 
5390 
$3) 
532. 
$33 
$36. 
$35! 
$35 
337 
$33 
639 
549 
$4} 
542 
543 
34 4. 
845: 
545. 
547 
#A9549 
#A%S49 
#A%559 
®#Aa%$5} 
#A®5S52 
#48553 
#A8S5564 
#A#S55' 
#a%5S56 
557 
$53. 
5593 
560 
$5) 
552. 
$$3 
564. 
$635! 
$65. 
$67 
$$8 
9593 
#579 
571 
$72 
$73 
574 
$735) 
575. 
577 
$73 
67>) 
989 
59} 
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580 
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640 


650 
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660 
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FORMAT (6F10,.0) 

FORMAT (6110) 

FORMAT (6F10.2) 

FORMAT (2F10-0) 

FORMAT (#NO OIL SANK WILL FORM BECAUSE THE. INITIAL! OLL SATURATION 
1 IS NOT GREATER THAN THE RESIOUAL OIL SATURATION. *) 

FORMAT (12F5.0) 

FORMAT (141) 

FORMAT (@6Xe*INTERNAL WELL BALANCE DEFEATED) 

FORMAY (8Xs#BOUNODING COMPLETE®) 

FORMAT (8Xe#RUNGACK COMPLETE®*) 

FORMAT (//o5Xe#WELL, NO 9225 AT #9E10e3e# o8E10039%.%) 

FORMAT (/910Xs*GAS RATE & #sE10e39%e OLL RATE = %yFB,28%- WATER R 
LATE = #e9FB,2,5%,%) 

FORMAT (/810Xs#CUM OIL = #9E10.39%0 CUM WAT @ %¢E10.38% GOR = %% 
1 El0e39% -o WOR ©& #5E10.3) 

FORMAT (/ISX,*FIELD: OIL RATE = @sE1Oed3e%s WATER RATE = #8£10,3) 
FORMAT (/I5Xe*FIELD CUM OIL 3 %6£10,39%s FIELO CUM WATER & @, 


1£10-¢3) 

FORMAT (/15Xe*TIME = HeELQe39*DAYS. CONDUCTIVITY RATIO = *9E£10,3 
1 ®e%s///) 

END 


SUBROUTINE TO BOUND IRREGULAR RESERVOIRS 
REF: JeeKe LINe PiHteDe DISSERTATIONs THE UNIVERSITY OF 
TEXAS AT AUSTIN (OECEMBERs1972) : 


SUBROUTINE S80UND 

COMMON /AAA/ Xw(125) 9 YWI125) oH (125)50(125) sx8(125) 5 
T¥B8 (125) sNEPsNPoNI s VOWet 

COMMON/BBB/ C(55255) 2D(55} 901 (55) 
DIMENSION XI (55) VI (55) 6XS8(1005 sySB(100) sA13(100) » 
LYI8 (i900) sAA(100955) 588(100) s9S(100) 
NBPI=NBP eo] 

PAI=£3,14]59 

LaNPeony 

00 lo J=lsN6u 

XJ] (Jp eXW(JeL) 

YIC(J) xv (Jeol) 

CONTINUE 

XB(NSP3) 2X84) 

YB (NBP1I)=¥S(1) 

CC20.5 

DO 29 I=IsN6P 

XM20,58(X8(Io1) oX8(7)) 

YMSO0 -S#(YB(Tel)*va(J)) 

DS(J=SORT (C(XB (Lol XB (7)) #829 (vYR( Tok) -¥B(1))*#2) 
HDS=cc#0s(T) 

TNEYB (Tel) oV3(T) 

TOFXB (141) X(T} 

THETASATANZ(TNs TO} 

BETAzTHETA<0,.S*PAI 

XIB(C1) eXneHDOS*COS (BETA) 

XSB(I) SAM*CHDOS*COS (BETACPAT) 

VYIB(Ij) eYRonOS*SIN( SETA) 

YS8(1) 2 YMeHOS*SIN(BETACPA) 

CONTINUE 

NaNGw 

MaNoP 

00 30 I71sN6P 

CO 30 J2IsNOw 
FNUMs(XSB(TI eX1( J) ) #82¢(YSE (EL) ~¥{ (J) ) 92 
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30 


40 


50 


60 
70 


80 


90 


Qa0aQ0 
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20 


30 


40 


FOENZ(XIB(I) eXI (J) ) @#2e(VYIB(T) VI (J)) *#2 
FACT=FNUM/FDOEN 

AA(19J) SALOG(FACT) 

CONTINUE 

DO 50 J=leNBP 

SUM=0.0 

DO 40 K=lolk 

FNE(XT3 (J) @Xw(K)) #26 (VIB (I) VW (Ky) #e2 
FD=(XS3 (J) XW (K) ) #820 (VSB) -YW(K) ) #42 
DSUM=sQ(K)/HI(K) FALOG(FN/FO) 
SUM=SUM*DSUM 

CONTINUE 

BB(J)=SUM 

CONTINUE 

DO 70 I=1eN 

DO 70 K=lel 

C(Kel)=0.0 

DO 69 J219M 

C(Ke IT) eC (Kol) *AA(UeT) MAAC JeK) 
CONTINUE 

C(loKk) eC (Kol) 

CONTINUE 

DO 8bO I=loeN 

0(1)=9.0 

DO 80 J=1eM 

O(T) =0¢(1) dB (U0) FAA UI) 
CONTINUE 

CALL SOLVE(N) 

DO 90 J=leNdw 

Q(J*L) SQ] (J) @H 

CONTINUE 

RETURN 

END 


SUBROUTINE SOLVE USED By SUBROUTINE BOUND » SE=: 30UND 
FOR REFERENCE. 


SUBROUTINE SDLVE(N) 
COMMON/BBB/ C€(55055) 00(55) #QI (55) 
DIMENSION A(55255) sID(55) 
NN=Ne} 
M2aN 
DO 10 I=loM 
A(IeNN) &0(T) 
IO0(1)<} 
DO 10 JEleN 
A(IeJ)=C(IoJ) 
CONTINUE 
Kz2] 
CONTINUE 
NR=K 
NC 8K 
BaABS(A(KyK)) 
DO 40 I=KoeM 
DO 40 JaKoeN 
IF (ABS(A(I9J))°B) 409940030 
CONTINUE 
NRsI 
NCaJ 
BaABS(A(IeJ)) 
CONTINUE 
IF (NRK) 70970950 
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50 


60 
70 


80 


90 


100 
110 


120 


1390 


150 


160 
170 
180 
190 
200 


CONTINUE 
DO 60 J2KsNN 
CmA(NRoVv) 
A(NRe J) =A(KeJ) 
A(KeJ)y eC 
CONTINUE 
CONTINUE 
IF (NC#K) 1009100980 
CONTINUE 
DO 90 I=l}9M 
CaA(I,NC) 
A(IleNC)=A(I eK) 
A(IeK)=C 
CONTINUE 
IsIO(NC) 
ID(NC) =IO(K) 
ID(K) se] 
CONTINUE 
CONTINUE 
IF (A(KeK)) 12091809120 
CONTINUE 
KK=Ke} 
DO 149 J=KKeNN 
A(KeJ) =A( Ke Jd) SACK 9X) 
DO 149 IT=loeM 
IF (Kel) 130912405130 
CONTINJE 
A(IeJ) HA(T oJ) eACI 9 K) PAC Key) 
CONTINUE 
K=KK 
IF (KeN) 2091109150 
CONTINJE 
DO 17y TeleN 
DO 170 J=leN 
IF (I1ID(J)e1I) 173991609176 
CONTINJE 
QI¢7)sACUeNN) 
CONTINUE 
Go TO 3190 
CONTINUE 
PRINT 200 
CONTINVE 
RETURN 
FORMAT (19H NO UNIQUE SOLUTION ) 
END 


THIS SUBROUTINE: PLOTS THE LOCATIONS OF THE OIL AND 
WATER FRONTS AS THEY DEVELOPE, ALL DATA IS TAKEN FROM 
THE MAIN PROGRAM. 


SUBROUTINE OBPLT 
COMMON/AAD/S XW(325) oe YW(125) 041 (125) 9Q(125) 9 XB(325) » 


1Y¥8(125) sNSPeNPoNIo NBwer 


COMMON/CCC/ X(80920) s¥ (80929) 9X0(80920) sYO(390920) » 


1WQ (80.20) sSCLeXMXeYMXeNSTeUMNP 9 RGO(80920) » 
CSGORsNOI We TMEI 


CALL SYMBOL (deS5e1lOeSeelSsT7HTIME = s27ue0s7) 
CALL NUMBER (999,.69999e09elG4e IMEI ,2o700080) 
CALL SYMBOL (999.09999209.14995H DAYS9270e0935) 
CALL PLT (Le l25*XMX#SCL 91,1257 YMXeSCL 9-3) 
CALL PLT (XB(NBP) ®SCL»sYK(NGOP) PSCL,3) 

00 10 I2leNsP 
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apt 
@pe 
ape 
ope 
ep. 
ope 
# pe 
% pt 
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ape 
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apd 
ape 
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20 
30 


40 


50 


CALL PLT (XB(1) *SCLeYB(1) #SCLe2) 
CONTINUE ae 
DO 30 Il» 
Q(I)-GE.0.0) GO TO 20 ; 
Pree cuse a atiescl + Vw cimeSCLs o7r1 903 1) 
GO TO 30 
E a 
falc ovaeen (XWCT) ®SCLeYWCT) FSCLe,0797909—1) 
CONTINUE a. 
20 JEleNPI 
4a IF (Q(J) 256.000) GO TO 50 
NSLEeNST#O (J) ZQMNP*0.5 
K2) eNSL ; 
Ee IF (ABS (XO0(J9K)) -GTeXMX) SO TO oh 
IF Pon ve econ oo P ie 
RGO(JeK) eLEe«SGO 
Pen ctor (XO 69K) SCL 9 YO(J9K) ®SCL 9007931908 1 
NTINVE 
i iF (ABS (X(JoK)) e¢GT.XMX) GO TO 50 
IF RIES crete co te Su 50 
WQ(J9K) -GT20-0) : 
Bi svtaol (XC Soh) OSCL oY (Ue K) 4SCp 90079409091) 
CONTINUE 
CALL PLT (02090-09999) 
RETURN 
ENO 


ope 
ope 
# pt 
ope 
opt 
#pe 
ape 
ape 
eps 
ope 
ope 
opt 
ope 
#Pe 
ope 
ape 
ope 
ope 
opt 


ope 


ope 
ope 
ape 
ope 
ape 
opt 
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13 
13 
20 
21 
22 
23 
24 
25, 
25 
27 
23 
29 
30 
31 
32 
33 
34 
35 
35 
37 
33 
39 
40 
41 
4a 
43 





elT75 


48, 
74 
500. 
710, 
"300. 
“1006 
90. 
"9206 
7006 
+5206 
3000. 
2500. 
2000¢ 
1500, 
1000. 
500, 
020 
“500° 
#10006 
“1500. 
-2000¢ 
~2300. 
25006 
~2800, 
©290V6 
-32006 
30006 
-3200. 
#3000, 
~3200, 
-3000. 
#3200, 
23000. 
#2750. 
2500, 
32150, 
-1800, 
1400, 
21000. 
-500, 
0.0 
500, 
1000. 
1600. 
20006 
25006 
30006 
3200. 
34006 
32006 
2800. 
26506¢ 
25006 
2250. 
20006 
1750, 
1500. 
1250. 


4 
6 


03S 


30, 
60006 
3006 
~340¢ 
T7006 
100. 
“5606 
4906 
“206 
-760¢e 
1500¢ 
17006 
1500. 
1700. 
1500¢ 
1700. 
15006¢ 
1700¢ 
15006 
1700-6 
1500¢ 
15506 
1300¢ 
13506 
1100¢ 
900e¢ 
600-6 
3006 
0,0 
5006 
1000, 
«1500, 
-2000, 
=2250, 
«2500, 
2600, 
ec€500, 
2600, 
=-2300, 
2100, 
“1700, 
#1350, 
-1000, 
-700¢ 
—300¢e 
100-¢ 
5006 
800, 
1100¢ 
13006 
1000¢ 
1000¢ 
10006 
10006 
1V0Q00¢ 
1000. 
10006 
1000.6 


40 


020 


500. 
S0.e 
“506 
41. 
“50 6 


~50 
4Q°« 
464 
3G0 
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4000. 
e75 
3e 

3e 

3e 


3e 


46, 
2.9 
1000-6 
1000. 
1000, 


1000. 
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1000. 10006 
750, 10006 
5006 1000¢ 
250, 10006 
0.0 10006 
“250. 1000¢ 
#5006 10006¢ 
750.6 10006 
©1000. 1000e 
“1250, 1000-6 
-1500. 10006 
“1/50. 10006 
“20006 1000¢ 
#2250, 1000. 
“2400-6 WV0e 
"25006 7506 
=2500+6 S00-« 
-2500, 250, 
"25006 0e0 
“25006 =“S00e 
“2500. ~-750e 
-2500, -l00U, 
72500. 1250, 
-2500, -1500, 
“2500-6 -1750, 
“2400. e1°00, 
“22506 “200U. 
“2000, 2000, 
~1750-¢ =2000, 
-1500, -c000, 
“1250, +1950, 
“10006 “1900. 
-7506 -1700,. 
500. -1500, 
=2506 #1350, 
0,0 1200, 
250.6 “1000, 
500, -B00e 
7506 “650¢ 
10006 -20Ce 
1250-6 =300e 
15006 “1006 
1750. 100. 
2000¢ 3006 
2150. 400, 
2300, 500, 
2500. 650, 
27006 BU0. 
2750, 900. 
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S0e 42.6 48. 

42< %4, 49. 

46. 46, 49, 

50 50 e 50. 

50, 50. 50. 

50, 50, SO. 

290071000000 1000.2000,3000, 
-10000,0 1000, 
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